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ABSTRACT 

Giant protoplanets evacuate a gap in their host protoplanetary disc, which gas must 
cross before it can be accreted. A magnetic field is likely carried into the gap, po¬ 
tentially influencing the flow. Gap crossing has been simulated with varying degrees 
of attention to field evolution (pure hydrodynamical, ideal, and resistive MHD), but 
as yet there has been no detailed assessment of the role of the field accounting for 
all three key non-ideal MHD effects: Ohmic resistivity, ambipolar diffusion, and Hall 
drift. 

We present a detailed investigation of gap magnetic field structure as determined 
by non-ideal effects. We assess susceptibility to turbulence induced by the magnetoro- 
tational instability, and angular momentum loss from large-scale fields. 

As full non-ideal simulations are computationally expensive, we take an a posteri¬ 
ori approach, estimating MHD quan tities from the pure hydrodynamical gap crossing 
simulation by Tanigawa et al. (20121. We calculate the ionisation fraction and estimate 
field strength and geometry to determine the strength of non-ideal effects. 

We find that the protoplanetary disc field would be easily drawn into the gap and 
circumplanetary disc. Hall drift dominates, so that much of the gap is conditionally 
MRI unstable depending on the alignment of the field and disc rotation axes. Field 
alignment also influences the strong toroidal field component permeating the gap. 
Large-scale magnetic forces are small in the circumplanetary disc, indicating they 
cannot drive accretion there. However, turbulence will be key during satellite growth 
as it affects critical disc features, such as the location of the ice line. 
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1 INTRODUCTION 

Giant planets capture their massive atmospheres from the 
surrounding protoplanetary disc. As the planet grows, its 
gravitational sphere of influence expands so that gas is cap¬ 
tured faster than the nebula can resupply it, and a gap 


is evacuated in the nebula around the planet (Lin & Pa- 


paloizou]|1985[ [Artymowicz fc Lubow|[l996| [Bryden et al. 
19991. Gaps have been seen in ~ /rm scattered light around 


TW Hya and infrared polarimetry of surface dust around 
HD97048 ( |Debes et al.|[?013| |Garufi et al.|[?014[ ). Spiral 
waves, likely launched by a protoplanet, have also been seen 
in HCO''" line emission within a gap in the HD 142527 cir- 
cumstellar disc ( [Casassus et al.|2013[ ). 

A circumplanetary disc encircles the planet following 


* E-mail: 

mark.wardle@mq.edu.au 


sarah.keith3@students.mq.edu.au; 


the collapse and detachment of the protoplanet envelope 
from the nebula ( Lunine fc Stevenson|[l982 Ayliffe & Bate 
2009a I. Gas flowing across the gap supplies the circum¬ 


planetary disc and controls the protoplanetary accretion 
rate (Lubow fc D’Angelo||200'6 1 . Gas with too much angu¬ 


lar momentum to reach the planet directly is captured by 
the circumplanetary disc and its flow is controlled by an¬ 
gular momentum transport processes operating within the 


circumplanetary disc (Fujii et al. 20111 |Lubow & Martin 

|2012t IRivier et al.| 

2012| 1 Tanigawa et al. 

2^^ hereafter 

TOM12; |Fujii et al. 

2014| [Keith & Wardle 

2^^ hereafter 

KW14; Szulagyi et al. 2014 Turner et al. 

20141. The cir- 


cumplanetary flow pattern potentially comprises of high- 
altitude inflow, and midplane Keplerian and outflow com¬ 
ponents (Bate et al. 2003 D’Angelo et al. 2003 TOM12, 


Gressel et al. 20131. Although circumplanetary discs are 


yet to be observed, there are prospects for their detection 
through infrared spectral energy distributions observed with 
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the James Webb Space Telescope or the Atacama Large 
Millimeter/submillimeter Array (ALMA; Isella et aL]|2014 


Dunhill|2015 Zhu|2015 1. This is strengthened by the recent 


observation of discs around three brown dwarfs, in infrared 
continuum and CO line emission ( Ricci et al.|2014 |. 

Magnetic helds likely play a role in the dynamics of 


the gap and circumplanetary disc system (Gressel et al. 


2013 Uribe et al. 2013 KW14, Turner et al. 20141. The 


field can have numerous effects, including generating hydro- 
magnetic turbulence via the magnetorotational instability 
(MRI; Balbus & Hawley 1991 Hawley et al.||19^ l, centrifu 


Wardle fc Konigl|1993 1, and magnetic braking (Matsumoto 


fc Tomisaka||200^ 


gaily driven disc winds and jets (Blandford &: Payne||1982 


Magnetic forces are transmitted by charged particles, 
and so these mechanisms require a minimum ionisation frac¬ 
tion. Collisions between charged and neutral particles give 
rise to non-ideal effects which counteract flux freezing - 
Ohmic resistivity. Hall drift, and ambipolar diffusion. Suf- 
hciently strong resistivity and ambipolar diffusion decouple 
the motion of the gas and field (see Turner et al.||20l4 and 
references within), and Hall drift can lead to complex held 


evolution sensitive to the global held orientation (Wardle 


19991 [Wardle & Salmeron||2012 

[Kunz & Lesur 

|20131 |Lesur 

et al.|2014[ O’Keeffe & Downes 

2014 Bai[[2015 



Gap-crossing dynamics is complex and so studies have 
relied on numerical simulations to model the how. These in¬ 


clude both grid-based (Bryden et al. 1999 Lubow et al.|1999 
Bate et al. 2003 Lubow & D’Angelo 2006[ TOM12) and 


smoothed particle ( [Bryden et al.|1999[|Ayliffe fc Bate|2009b[ 


20121 hydrodynamical simulations, and ideal magnetohydro- 
dynamical (MHD) simulations (Nelson fc Papaloizou|[2003 


jPapaloizou et aL]|2004| [Uribe et al.||2013| ). Global 3D MHD 

simulations of the wider system have also been used to gen¬ 


erate synthetic observation maps for ALMA (Flock et al. 

[Ms] ). 

The recent inclusion of Ohmic resistivity adds an im¬ 
portant level of realism to gap-crossing modelling ([Gressel[ 


et al. 20131. Ohmic resistivity is strong in the dense cir¬ 


cumplanetary disc, and can produce an MRI stable, dead 
zone with implications for planet growth and disc evolution. 
However, non-ideal effects are also at work in relatively low 
density regions where the Hall effect and ambipolar diffusion 
can be strong ( [Wardle|[2007[ ). Their inclusion may enhance 
or suppress turbulence, affecting the how and protoplanet 
growth rate.Indeed, the addition of ambipolar dihusion to 
resistive simulations radically alters the protoplanetary disc 
how between 1-5 au, producing laminar accretion powered 
by a magnetocentrifugal wind, rather than the traditional 
accreting turbulent surface layers ( [Gressel et al.[[^15[ ). 

In this article we perform a comprehensive study of the 
relative importance of the three non-ideal effects on huid dy¬ 
namics throughout the gap. Owing to the signihcant com¬ 
putational cost of including all three ehects in a full, 3D 
simulation, we take an a posteriori, semi-analytic approach. 
We calculate MHD quantities using a snapshot from the 
TOM12 3D hydrodynamical simulation (described in Sec¬ 
tion EH- This allows us to calculate detailed ionisation 
maps including ionisation from cosmic-rays, stellar X-rays 
and radioactive decay, accounting for grain charging (Sec¬ 
tions [2.2[ and ihl We determine the strength of non-ideal 
effects (Sections 2.3 and 3.21 to estimate the held strength 


and geometry (Sections 2.4 and 3.31. We hnd that the held 
would not alter the general huid motion, justifying our use 
of an underlying hydrodynamical gap-crossing model (i.e., 
TOM12). The implications of incorporating magnetic helds 
and non-ideal effects in gap crossing, such as the additional 
force from large-scale helds, are considered (Sections[4.1[and 


4.2 I. Finally, we present a summary and discussion of hnd- 
ings in Section]^ 


2 MODEL DESCRIPTION 

In this section we outline the disc model used to describe 
the protoplanetary disc and gap. We give details of the ion¬ 
isation and dihusivity calculations, along with estimates for 
the magnetic held strength and geometry. 


2.1 Disc and gap structure 

We take a semi-analytic, a posteriori approach to the calcu¬ 
lation. We use a snapshot of a pure hydrodynamical simu¬ 
lation as the basis for estimating MHD quantities. The pro¬ 
toplanetary disc and gap model that we use is the TOM12 
three-dimensional, hydrodynamical simulation. The simula¬ 
tion models a protoplanetary disc surrounding a star of mass 
M*, in which a gap has been carved out through gas cap¬ 
ture by an embedded protoplanet of mass Mp, at hxed or¬ 
bital radius dp. The simulation corresponds to a protoplanet 
with the present-day mass and orbital radius of Jupiter (i.e., 
Mp = Mj, and dp = dj, where dj = 5.2au). Gas is treated 
as inviscid, and self-gravitational and magnetic forces are 
neglected. 

The planet is located at the origin of a local cartesian 
co-ordinate system {x,y,z), which orbits the star at an or¬ 
bital angular frequency Qp. The x-axis is oriented along the 
radial direction, r, which extends from the star to the planet; 
the j/-axis is oriented in the azimuthal direction, cf), parallel 
to the the protoplanet orbit; and the z-axis is parallel to 
the protoplanetary disc angular momentum vector. To cap¬ 
ture details of the Hne flow structure near the protoplanet 
the simulation has eleven levels of nested grids, each with 
rij; X Uy X — 64: X 64: X 16 data points. The total simu¬ 
lated portion of the disc extends over x G [—12Hp, 12Hp], 
y G [—12Hp/dp, 12Hp/dp\, and 2 G [0, 6LIp] where Hp is the 
scale height of the protoplanetary disc. The local cartesian 
approximation requires that the Hill radius, 

Rh = 


which is roughly the feeding zone of the protoplanet, satisfies 
Rh <C dp. Here Mq is the mass of the Sun. The simulation 
data used here was taken after 160.7 protoplanet orbits, al¬ 
lowing the simulation to approach a steady state. The gap 
reached a maximum density contrast relative to the unper¬ 
turbed disc of ~ 3. 

The simulation used a non-dimensionalised form for 
the MHD equations and so we must rescale the velocity 
V = (vj;,Vy,Vz) and density p data for our calculations. 


M 


(sM* 


0.36 au 


d \ / M y 

5.2au ) \Mj) 


M, 


( 1 ) 
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We rescale the column density using the minimum mass So¬ 
lar nebula column density at the orbital radius of Jupiter 
So = 


140 gem ^ iKitamura et al.]|2002 Williams & Cieza 


20111. This model estimates the gas mass in the Solar neb¬ 


ula by augmenting the solid material in Solar System planets 
with sufficient hydrogen and helium to bring the composi¬ 
tion up to Solar | Weidenschilling||1977 Hay^hi||1981 Ob¬ 
servations of Class II Young Stellar Objects (i.e., those with 
discs and strong UV and Ha emission indicating active ac¬ 
cretion) are currently limited to ~ 20 au resolution, however 
inferred disc column density profiles are broadly consistent 
with this profile. As the gap gas density naturally drops over 
time, we allow for simple reduction of the column density, 
E(a;, y, z), though multiplication of a constant parameter /s, 
where /e = 1 recovers the original TOM12 results. 

The simulation is isothermal, and we adopt the tem¬ 
perature of a black body in radiative equilibrium with So¬ 
lar luminosity, excluding any heating from inflow processes, 
at the orbital radius of Jupiter, T = 120 K ( Wardle||2007 |. 
This is supported by interferometric CO images of T Tauri 


and dust temperatures from SED modelling (Guilloteau & 


Dutreyj 1998 Andrews fc Williams|2005 Pietu et al.|2007 |. 

We adopt Solar gas composition, with a mixture of 
80% molecular hydrogen and 20% atomic helium. This cor¬ 
responds to a mean molecular weight of = 2.34mp, 
where rup is the mass of a proton. This allows us to calcu¬ 
late the neutral number density, n, isothermal sound speed 
Cs = 0.65kms“^, and aspect ratio: 


— = — = 4 5 
d Vk 


X 10"^ — 


( 2 ) 


where Vk = (GMt/d)^ is the Keplerian velocity. 

The flow pattern within a gap transitions from orbit¬ 
ing the star [orbital radius (x + dp)^ -|- y^] , to orbiting the 
planet (orbital radius yjx'^ -\- y^) within the Hill sphere. In 
our analytic calculations we approximate the flow geometry 
by a composite Keplerian angular velocity functions: 


Q 


^ f [cMp {x^ + y^) 3] ' for + y^ < R„, 


[^GAf+(x J- dp 


(3) 


otherwise. 


TOM12 also simulated the circumplanetary disc, how¬ 
ever they note that as it is shear-dominated, artificial vis¬ 
cosity is strong and so the results are more reliable in the 
gap. Similarly we include the circumplanetary disc in our 
calculations with the caveat that the disc structure in this 
region is uncertain. 

The scale height in the circumplanetary disc. 


H, = Hp 


\/x'^ -b y'^ 

3Rh 


5.2 X 


\/x‘^ + y^ 
0.2 Rh 


(4) 


(5) 


is significantly lower, and this needs to be accounted for 


when calculating gradients in the fluid flow (see Section 2.31. 
To that end we use a composite scale-height, H, which tran¬ 
sitions at the boundary of the circumplanetary disc hydro¬ 


static region: 

He for < 0.2Rh, z/ < 5, and 


H = 


Hp 


otherwise. 


( 6 ) 


2.2 Degree of ionisation 

Charged particles transmit magnetic forces to the bulk neu¬ 
tral flow. In the gap, ionisation is largely caused by penetrat¬ 
ing cosmic-rays and stellar X-rays, however there is also a 
weak, pervasive contribution from decaying nuclides. We cal¬ 
culate the ionisation level by solving the coupled rate equa¬ 
tions for electron and ion number densities and mean grain 
charge subject to overall charge neutrality. We follow the 
method outlined in KW14, but using a more detailed cal¬ 
culation of the grain charge. A summary of the method is 
given below. 

We calculate the grain number density, Ug , for spherical 
grains with radius Ug = 0.1/rm and bulk density 3gcm“® 
with a constant dust to gas mass ratio fdg ~ to ac¬ 

count for incorporation of solids into protoplanetary bodies 


(Pollack et al. 19941. 


Internal ionisation from radioactive decay is primarily 
from the short-lived radionuclide ^®A1, which c ontributes an 
ionisation rate Cr = 7.6 x 10“^® s“^ {fdg/ fo) (Umebayashi 
& Nakano 2009| ). The ionisation rate has been scaled by 
the ratio of the dust mass fraction to the Solar photospheric 
metallicity, /q = 1.3%, so that the ^®A1 abundance is consis¬ 
tent ( Asplund et al.|2009 |. The midplane dust fraction will 
be enhanced by settling, however the consequent increase in 
the ionisation fraction will be small as ionisation is predom¬ 
inantly by external radiation. 

External ionisation sources are attenuated by the col¬ 
umn density above and below each location in the simula¬ 
tion. Cosmic-ray ionisation occurs at the interstellar cosmic- 
ray ionisation rate, ^cr = 10“^^ s“^, attenuated by the at¬ 
tenuation depth Ecr = 96gcm“^ approaching from above 


and below the disc (Umebayashi & Nakano 1981 20091. 


We account for ionisation from diffuse scattered stellar X- 
rays but neglect direct illumination, instead assuming that 
shielding by the inner portion of the disc is effective. Ion¬ 
isation for solar luminosity stellar X-rays occurs at a rate 
CxR = 9.6 X 10~ ^^ s~^ (d/dj)~^, with attenuation depth 
Exr = 8gcm“^ (Igea & Classgold 


1999 


Turner & Sano 


20081. The total ionisation rate, i(, is the sum from the three 


sources C = Cr + Ccr + Cxr- 

We solve for the electron, ion and grain charge number 
densities by calculating the grain charge needed for overall 
charge neutrality in the steady state [see equations (26)- 
(29) in KW14[. We include charge focussing in the electron 
and ion capture rate coefficients for grains (Umebayashi & 
Nakano||1980 1: 


kig - 'XQjg 


8kbT j Zg<0 


'KVCli 


, 2 /SfciT 

keg = Tvagyl - X 


^ agk^T ) 

^ ^ agktT J 


Zg>0 


Zg < 0 


1 exp(-f^) ^g>0 


(7) 


( 8 ) 


where Zgq is the mean grain charge, kb is Boltzmann’s con¬ 
stant, and me is the electron mass. We model metals as a 
single species, adopting the mass, rrn, and abundance, Xi, of 
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the most abundant metal - magnesium ( Lide|20d4 Asplund 
et al. 12009 \. The recombination rate for ions and electrons is 

* i<-». _ n'7 o 1 i_ __ 


fcei = 1.20 X 10 (T/IOOOK) cm^s ^ (McElroy et al. 


20131. If the grain charge is known these equations give the 


ion and electron number densities: 


‘2kei 

ke 


1 + 


4:kein( 

kegkigTl^ 


- 1 


Tie 


kig 


(9) 

( 10 ) 


Equations 0 -( |10[ ) are analogous to equations (30)- 
(33) in KW14 except for the inclusion of a non-zero grain 
charge for charge focussing. The protoplanetary disc has 
a higher average ionisation fraction resulting in a higher 
grain charge which cannot be neglected (as it could be 
in the KW14 circumplanetary disc analysis). We deter¬ 
mine the equilibrium values of ne,ni, and Zg numerically 
by solving equations ([7|-( |10[ ) along with charge neutral¬ 
ity, for the grain charge. We use the Brent-Dekker method, 
gsl_root_fsolver_brent, for root-finding implemented in 
the GNU Scientific Library to solve for the grain charge to 
an accuracy of 1% ( [Galassi et al.|[^09[ ). This method com¬ 
bines bisection and secant methods for rapid convergence. 

We also determine the effectiveness of collisional ionisa¬ 
tion produced in MRI turbulent regions. Gurrents generated 
by an MRI held may be accelerated sufhciently to ionise neu¬ 
tral particles, providing additional ionisation ||Muranushi| 


et al. 2012 Okuzumi & Inutsuka 20151. We calculate the 


kinetic energy of accelerated particles using equations (37)- 
(38) in KW14 to determine if it exceeds the ionisation po¬ 
tential of neutral species (e.g., potassium has the lowest ion¬ 
isation energy of the abundant elements, with E ~ 4.34 eV). 
We hnd that although electrons reach energies of up to 1 eV 
in the disc atmosphere {z ~ 5Hp) where collisions are infre¬ 
quent and the mean-free path is long, in the denser midplane 
electrons are not accelerated above 0.05 eV. This is too low 
to ionise any atomic species and so the MRI does not con¬ 
tribute signihcantly to the ionisation. Thus, we neglect the 
effect in our calculation of the ionisation fraction. 

Finally, we note that thermal ionisation is likely not ef¬ 
fective in the gap surrounding Jupiter. The protoplanet ary 
disc at d — dj is much cooler than the required T ~ 1000 K, 
and further shock heating of accreting material at the cir¬ 
cumplanetary disc is weak. In the circumplanetary shock, 
shock heating flux Fb ~ a{40K)‘^ (calculated for an accre¬ 
tion rate 10“®Mj/year, at r ~ 300-Rj from the protoplanet; 
[Cassen fc Moosman|198l| is negligible compared with Solar 
irradiation, = ct( 120K)‘^ at this orbital radius, preclud¬ 
ing thermal ionisation. 


where rjo, Vh, and rjA are the transport coefficients for 
Ohmic resistivity. Hall drift, and ambipolar diffusion. These 
non-ideal effects cause the magnetic field to slip through the 
gas, and can dominate the field evolution if they are too 
strong. They are caused by neutral collisions disrupting the 
E X B drift of charged particles. Differences in the density 
dependence of these effects leads to three different regimes 
in which Ohmic resistivity. Hall drift and ambipolar diffu¬ 
sion are most effective in high, moderate, and low density 
regions, respectively. For a predominantly vertical magnetic 
held, Ohmic resistivity and Ambipolar diffusion act similarly 
(and distinctly from the Hall drift), and so they are often 
combined into the Pedersen diffusion rjp = rjo +riA- We cal¬ 
culate the strength of the three effects using the transport 
coefficients from equations (53)-(55) in KWl-^ 

How strong these effects must be to dominate the in¬ 
ductive term and decouple the motion of the held and gas 
depends on the held geometry. Below we outline the cou¬ 
pling conditions for MRI turbulent, toroidal, and poloidal 
components in turn. 

The MRI dynamo harnesses shear to generate a tur¬ 
bulent held. Dihusion restricts MRI turbulence by limiting 
bending of the smallest-scale held modes. The MRI oper¬ 
ates effectively if the fastest growing mode, with wavelength 
A = 2TvVa/^, survives diffusion. Magnetic diffusion counter¬ 
acts dynamo growth of short wavelength huctuations with 
wavelength A < rj/ra,z, for which the diffusion rate exceeds 
the growth rate. Here Va,z = Bz/yJ^-np the vertical Alfven 
speed. Thus, the MRI operates in effectively ideal MHD con¬ 
ditions if Pedersen and Hall components are simultaneously 
below the threshold (Sano fc Stone|2002 Turner et al.|2014|: 


Vp ^ ^ and \r]H 




( 12 ) 


Here, the absolute value of the Hall term rehects its signed 
nature. The MRI also requires the held must be weak enough 
that the fastest-growing mode is conhned within the disc 
scale height, so that for ideal MRI ([Okuzumi fc Hirose|201l[|: 


I3.= 


S-KpCg 


> 8^-^ 


(13) 


Unlike Ohmic resistivity and ambipolar diffusion. Hall 
drift is not diffusiv^ and can change the behaviour of MRI 
turbulence, even in the presence of strong Ohmic and am¬ 
bipolar dihusion. Hall drift is antiparallel to the current den¬ 
sity and may cooperate with, or act against Keplerian shear 
according to the relative orientation of the protoplanetary 
disc angular momentum vector and vertical component of 
the magnetic hel d, s = sign(B ■ fl) (jWardle 1999 Balbus 
fc Terquem||200f Wardle & Salmeron 20121. For example. 


2.3 Non-ideal MHD effects and magnetic coupling 

The evolution of the magnetic held and its role In the gas 
dynamics depends on the strength and nature of the mag¬ 
netic coupling between the held and gas. Magnetic coupling 
is manifest in the induction equation, 

— = V X (v X B) — V X [7?o(V X B) + 7?j/(V X B) X B] 


-V X [77a(V X B)_l], 


( 11 ) 


^ There is a typographical error in the expression for Ohmic 
resistivity in equation (53) of KW14. The scaling of the 
diffusivity with the neutral density should be r}o oc p 
rather than the inverse proportionality given. Additionally, the 
grain Hall Parameter in equation (51) of KW14 is a fac¬ 
tor of yjvriglmn too large, and should be written as 0g = 
5.5 X 10-^Zg {B/1 G) (10^^ cm-3/n) (0.1 p-m/agf >/l000 K/T, 
instead. This reduction has no impact on the KW14 results, as 
^g was al ready so small as to b e negligible. 

^ Indeed, |Kunz &: Lesur| | |2013[ | found that the Hall effect can be 
anti-diffusive in unstratified shearing boxes. 
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if the field is aligned with the rotation axis, the Hall ef¬ 
fect can destabilise the flow, enhancing turbulence. On the 
other hand, if the field and rotation axis are anti-aligned, 
Hall drift acts against the shear, suppressing turbulence. In 
the Hall-MRI regime, the Hall effect counteracts Pedersen 
stabilisation against field tangling if it dominates and the 


vertical field orientation is favourable (Wardle & Salmeron 
2012 Lesur et al.|2014 |. Hall-MRI requires 

\riH\ > Vp, \Vh\ > and s = 1, (14) 

along with a Hall-MRI weak-field criterion, which we show 
in Appendix|^is the same as equation (131. Indeed, the Hall 
effect need not dominate diffusion to impact turbulence, as 
even weak Hall drift modifies the wave number and growth 
rate of the fastest growing mode ( Wardle fc Salmeron|2012' I. 

A toroidal field surrounding the star or protoplanet, 
must also be maintained against diffusion which would tend 
to unwrap and straighten field lines. The impact of non¬ 
ideal effects on a toroidal filed are most important close to 
the central object where gradients are strongest. A toroidal 
field is preserved against diffusion if magnetic induction ex¬ 
ceeds non-ideal effects, as captured by the coupling threshold 
(c.f.. Turner fc Sano|2008| see Appendix 0 : 

(15) 


rtp < and |? 7 h| < 


A toroidal field can also be preserved against strong dif¬ 
fusion if Hall drift dominates. Hall drift is along the current, 
which is in the fluid rotation direction for both the aligned 
and anti-aligned cases (i.e., s = 1 and s = — 1). So, either by 
enhancing the field wrapping (i.e., for s = 1 ) or by unwind¬ 
ing and rewrapping the field in the counter flow direction 
(i.e., for s = —1), the Hall effect produces a strong negative 
toroidal field. See Appendixj^for a discussion of the effect of 
vertical field direction dependence for coupling of a toroidal 
field. Simulations are needed to verify this behaviour, but 
we highlight Hall-dominated regions with the anticipation it 
may extend the reach of the toroidal field. We simply require 
that azimuthal winding from Hall drift exceeds diffusion: 

IrjHl > Tjp, \r]H\ > (16) 


Finally, the disc contains a poloidal field component 
throughout. Although magnetic coupling limits the mini¬ 
mum field bending radius of the poloidal component, an es¬ 
sentially vertical field will always couple, owing to its small 
gradient. As we are not concerned with the exact geometry 
of the poloidal component, we simply assume that it is ap¬ 
proximately vertical, couples everywhere and so permeates 
the gap. 


2.4 Magnetic field strength and geometry 


Protoplanetary discs inherit a large-scale magnetic field from 
their progenitor molecular cloud. Field measurements are 
difficult as current observations are limited to ~ 100 au res¬ 
olution ( [Stephens et al.|2014[ ). The disc field is certainly com- 
pressionally enhanced over the cloud field {B ~ 1-100 mG; 


Shu et al. 19871, while the equipartition field Heq = Snp pro¬ 


vides a maximum field strength the disc will support before 
magnetic forces exceed the thermal pressure, p = pc^- 

Nevertheless, it is possible to gain a more precise esti¬ 
mate from the stellar accretion rate, M ~ lO“®M 0 yr“^, 
as the field is thought to play a principle role in angular 


momentum transport (Hartmann et al. 1998 see references 
McKee fc Ostriker|2007|. This can be made from the az¬ 


imuthal component of the axisymmetric momentum equa¬ 


tion, vertically integrated between the disc surfaces (Wardle 

(B ■ VB), 


20071, 


P[(v'V)v]^ = 


47r 


(17) 


This yields a minimum field strength for vertical, toroidal, 
and turbulent components in the disc ( Wardle|2007l: 


Hz = 


6.0 X 10"^ G 

n 

X 


M 


10-® Mq yr- 


Bs — 


1.7 X 10-8 s-i 
2.8 X 10“^ G 
Q 


M 


10-8 Mq yr- 


X 


1.7 X 10-8 s- 


Hr, 


0.23 au 


(18) 


and (19) 


Bmri = 


5.7 X 10"^ G 


n 


M 


10-8 Mq yr-1 


1.7 X 10-8 s-i 


Hr, 


0.23 au 


( 20 ) 


respectively. Components of the turbulent field are (Sano 
jet al.|2004D : 

Hmri.t = 0.35Hmri, (21) 

Hmri.c/) = 0.92Hmri, and ( 22 ) 

Hmri.z = 0.19Hmri. (23) 


We now turn to estimating the magnetic field strength 
in the gap. Turbulence would be continuously generated 
there so that equations ( | 20 | )-l | 2 ^ remain valid, however our 
estimates of the vertical and toroidal fields are not sufficient 
in the gap. 

The vertical component would be drawn in with the flow 
and reduced or enhanced in keeping with the large density 
variation across the gap. We assume that the field is frozen 
into the gas at the outer edge of the gap (at Xg = 1.45 au), 
drawn into the gap by the flow. The self-consistency of our 
adoption of this assumption is verified by our results (see 
Sections 3.2 4.2 1 . The field would be pinned to the flow at 


the lowest magnetically-coupled gas layer, reducing or en¬ 
hancing the field according to H oc E-^. Nevertheless, for 
simplicity, we assume the field is pinned at the midplane, 
with a field strength B^^g = 4.4 x 10“^ G, given by equation 
(18l at Xg = 1.45au, j/ = 0 where Eg = 97gcm-^. Normal¬ 
ising to the field strength at the outer edge of the gap, this 
leads to a vertical field strength profile of 


Hz = Hz 


(24) 


in the gap. 

Determining the toroidal field strength is considerably 
more difficult as the field would be wound up continually 
until field-line drift balances differential rotation. As this 
is beyond the scope of our treatment, we simply assume 
the constant ratio y/djHp between the toroidal and vertical 
components in the disc is preserved for our estimate in the 
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gap. This leads to a gap toroidal field strength proHle of: 


Btj, 





(25) 


where = 2.0 x 10~^ G is the toroidal field strength, 
likewise taken at Xg = 1.45 an, y — 0. We cap the held 
strength so that the associated magnetic pressure does not 
exceed the gas pressure: 


SyrpCs 


> 1 . 


(26) 


This allows the held to expand in response to strong mag¬ 
netic pressure at high altitude, correspondingly reducing the 
held strength. 

Numerical simulations suggest that the toroidal com¬ 
ponent may be subject to periodic polarity reversals, as 
it is lost through magnetic buoyancy every ~ 40 orbits 


(Miller & Stone 2000 Shi et al. 2010 Hirose & Turner 


2011 Flock et al. 20121. The magnitude of the toroidal 


held would also vary vertically between its surface value, 
B^^s, and zero at the midplane, with the approximate scal¬ 
ing B^ ~ B^^s[z/Hp), appropriate for a thin, rotationally- 
supported disc that is well coupled to the magnetic held 
( Wardle fc Koni^|1993 |. For simplicity we take con¬ 
stant with height, but account for vertical gradients (e.g., in 
Appendix]^, and probe the effect of uncertainty by global 
enhancement or reduction of the the poloidal and toroidal 
components through multiplication by a constant parameter 
/s- 


2 . 4.1 Self-consistent coupled field geometry 

The variation of diffusion and held gradients across the gap 
means that not all three held components are present ev¬ 
erywhere. Here we determine where diffusion permits each 
held component. We use the diffusivities to determine the 
magnetic structure based on the level of coupling for each 
held component (poloidal, toroidal, and turbulent), and only 
include the components where they couple. Our procedure 
for calculating the held geometry is: 


(i) Poloidal held - we include a poloidal component at all 
locations, calculated with equation ( |18| l in the main proto¬ 
planetary disc (i.e., |a;| > Xg), or using the hux-conserved 
form, equation (241, inside the gap. It is safe to include a 
poloidal component everywhere since a pure vertical held 
always couples and that any radial bending will adjust to 
the level of diffusion. 

(ii) Toroidal held - we calculate the toroidal component 
strength using equation (191 in the bulk protoplanetary disc 
how, and equation (251 in the gap. We cap the held strength 


by the gas pressure through equation ([26f. We calculate rjo, 
rjH, and pA using the net held strength, B = (B^ -I- B^)^^^, 
and determine if the held is coupled using equations (|15[) 


and (I 61 . If the held is coupled the toroidal component is 
kept; otherwise it is set to zero. 

(iii) MRI turbulent held - hnally, we determine if an 
MRI turbulent held is sustained. We calculate the tur¬ 
bulent held strength in both the protoplanetary disc 
and the gap using equations ( |20[ )-( [23| |. We calculate 
Vo, Vh, and tja using the net held strength, B = 


[-^MRi,r 4“ {Bfi, Bmri,(/))^ + i.Bz T Bmri.z)^] , and de¬ 
termine if the held is sustained using equations (12 1 -( 13 1 . 


3 RESULTS 

Here we present the results of the calculations developed in 
Section]^ including the ionisation fraction, magnetic held, 
and strength of non-ideal MHD effects. Figures are shown 
for a Jupiter-mass planet, orbiting a Solar mass star, at the 
present orbital distance of Jupiter. Unless otherwise stated, 
we take fdg = 10 “"^, /e = 1 , /b = 1 , and s = 1 . 


3.1 Degree of ionisation 

Fig.[T] shows the ionisation prohles within the gap and disc. 
The left and centre panels show 2D slices through the gap, 
with the colour scale showing the ionisation fraction. The 
left panel shows the top-down view in the x-y plane at the 
midplane, and the centre panel shows the edge-on prohle 
in the x-z plane at y — 0. The star is located beyond the 
simulation boundary at a: = —dj, y = 0. Small white patches 
in the top- and bottom-right corners of the edge-on prohles 
are beyond the simulation boundary. The right panel shows 
the dependence of the vertical ionisation prohle on the dust 
to gas mass ratio, at a; = 0.3 au, y — 0. 

Shielding from the overlying gas reduces the ionisation 
fraction with height and in denser regions. X-rays penetrate 
to 2 ~ Bp, while cosmic rays reach everywhere except for 
the circumplanetary disc (located at x,y < 0.1 au), owing to 
its high column density. Instead, ionisation in the circum¬ 
planetary disc is from radioactive decay, which is a much 
weaker ionising source. Consequently, the ionisation frac¬ 
tion in the circumplanetary disc is much weaker than in the 
protoplanetary disc, as supported by previous studies (|Fu-| 
jii et al.|[2MI| [Mil KW14, [TuFner et al.|[M4| ). The sharp 
drop in the ionisation fraction prohle at 2 ~ Bp occurs at 
the shock boundary surrounding the circumplanetary disc, 
visible as a halo of poorly ionised gas in the edge-on view 
shown in the centre panel of Fig. 

Dust grains reduce the ionisation fraction by soaking up 
free electrons. Grain charge capture is signihcant and most 
efficient at the midplane where grains are abundant. It also 
reduces radioactive decay, inhuencing ionisation, and hence 
the potential for magnetic coupling, in the circumplanetary 
disc. 


3.2 Non-ideal MHD effects and magnetic coupling 

Fig. [2] shows diffusivity prohles for the gap. The left and 
centre panels show 2D slices through the gap, with the 
colour scale showing the held perpendicular diffusivity r]± — 
{rip , along with logarithmically spaced contour lev¬ 

els. The left panel shows the top-down view at the midplane, 
and the centre panel shows the edge-on prohle at j/ = 0. The 
right panel shows the vertical prohles of the transport coef- 
hcients of Ohmic resistivity, the Hall effect and ambipolar 
diffusion, along with the coupling threshold for toroidal and 
MRI helds. 

Ohmic resistivity traces the high density structures such 
as the spiral arms and circumplanetary disc and is lowest 
in the evacuated regions to the upper-right and lower-left of 
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logio(ng/n) logio{ng/n) z [H] 


Figure 1. Ionisation fraction profile slices through the gap, a posteriori computed for the hydro dynamical simulation of |Tanigawa et ah] 
( |2012| |. Left and centre panels show 2D slices of the ionisation fraction logiQ(ne/'n) as contour plots with logarithmically spaced contours 
for y = 0 and 2 : = 0, respectively, and for f^g = 10““^. Right panel shows nejn as a function of height for x = 0.3au, y = 0 for dust to 
gas mass ratios f^g = 0,10““^, 10“^ with solid, dashed, and dotted curves respectively. 



i|ogio(tip +nH ) [cm^s-^ 2'°9io(np +% ) [cm^s ^ 



z[H] 


Figure 2. As for Figure but showing a contour plot of the field-perpendicular diffusivity, 77 ^ = ydp + ’Ih the midplane (left 
panel), and y = 0 (centre panel). Ohmic resistivity (solid curve), Hall effect (dashed curve), and ambipolar diffusion (dotted curve) are 
also shown as functions of height at x = 0.3 au, y = 0 (right panel). The long-dashed and dot-dashed curves show the coupling thresholds 
for toroidal and MRI fields, and respectively. 


the circumplanetary disc. It decreases with height in keeping 
with the exponential density profile. Hall drift is relatively 
constant with height above the midplane and dominates up 
to 2 ; = 2.5 Hp. Ambipolar diffusion is strongest in the low 
density atmosphere where neutrals are too sparse to influ¬ 
ence electron and ion motion. 

Non-ideal effects are orders of magnitude below the level 
at which the toroidal field decouples, except in the circum¬ 
planetary disc. Here the field is uncoupled to the midplane 
gas flow, but is anchored to, and transported by, gas in the 
overlying coupled region. 

Non-ideal effects are below the MRI suppression thresh¬ 
old in the disc atmosphere, 2 > Hp. The turbulent region 
is limited by the weak field condition [equation ( |13| l], rather 
than diffusion. Magnetic pressure is weaker than gas pres¬ 


sure at 2 < 2Hp, and the field is weak enough to bend. 
The strong field region reaches towards the midplane near 
the planet where field enhancement from flux-conservation 
is greatest. 


3.3 Magnetic field strength and geometry 

Fig. [3] shows the magnetic field strength and geometry. The 
top panels are contour plots, in the top-down (top-left panel) 
and edge-on (top-right panel) views, with colour scale show¬ 
ing log,Q(H), calculated for s = 1. The bottom-left panel 
shows vertical profiles for the poloidal, toroidal and turbu¬ 
lent magnetic field components evaluated at 2 ; = 0.3 au, y = 
0. A turbulent field is only present in the active zone be¬ 
tween 0.5 < 2 < 2 Hp for s = — 1 and for 2 < 2Hp for 
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Figure 3. Properties of the estimated magnetic field in the gap for the standard parameter set = 10“^, /e = 1, /s = 1. Top-left and 
-right panels show contour plots of the field strength, B, at the midplane and y = 0, respectively. Bottom-left panel shows the vertical 
profile of field components at r = 0.3 au, y = 0. The vertical, toroidal, MRI field with s = 1, and MRI field with s = — 1 are shown as the 
dot-dashed, dashed, solid, and dotted curves respectively. Bottom-right panel shows the edge-on view of the field geometry, colour-coded 
according to the dominant field component for s = 1: (a) poloidal field (yellow), (b) toroidal field (green), and (c) MRI turbulent field 
(blue). Orange, dark green, and dark blue regions indicate where the Hall effect sustains and influences the toroidal field. Hatched MRI 
unstable regions show where Hall effect dominates, so that the field is predominantly turbulent if s = 1, and toroidal if s = —1. 


s = 1. The bottom-right panel is an edge-on view of the gap, 
colour coded according to the dominant field component. 
The gap Is divided into the following regions: (a) poloidal 
(yellow), (b) toroidal (green), and (c) MRI turbulent (blue) 
field. Orange, dark green, and dark blue regions indicate the 
corresponding regions where the Hall effect maintains cou¬ 
pling of the toroidal component, and so the field may be 
counter-wrapped when is aligned with the rotation axis. 
Hatched blue regions are Hall-MRI unstable and so they are 
predominantly turbulent if s = 1, or predominantly toroidal 
If s = —1. 

The field strength is relatively uniform in the gap but 
flux conservation enhances the field in the circumplanetary 
disc. A toroidal field would easily couple and permeate the 
gap, but is limited by magnetic pressure in the disc atmo¬ 
sphere, and so the field would be increasingly poloidal with 
height. The gap would be turbulent between z ~ 0.5-2 Hp, 
except in the circumplanetary disc where the field is too 
strong. The Hall effect extends the turbulent region to the 
protoplanetary disc and circumplanetary disc midplanes if 
s = 1, however the midplanes are MRI stable if s = —1. The 
upper boundary of the turbulent zone is limited by the weak 
field requirement, rather than non-ideal effects. 


Fig. a shows the effect of varying parameters on the 
field geometry: fdg = 0, 10“^ (top- and bottom-left panels), 
/e = 0.1, 10 (top-and bottom-centre panels), and fs = 0.1, 
10 (top- and bottom-right panels). In general, we find that 
the gap is either predominantly toroidal, ideal MRI, or Hall- 
MRI unstable. The upper MRI boundary only depends on 
Pz through the column density and magnetic field strength, 
whereas the boundary of the Hall-MRI layer is controlled by 
the height at which the Hall drift exceeds the MRI threshold. 
Thus, the Hall-MRI region extends by lowering the 
ionisation fraction, caused either through enhanced surface 
density attenuating ionising radiation (e.g., for /e > 1 ), or 
by grains soaking up free electrons (e.g., for fdg = 10“^). 
Lowering the magnetic field strength reduces the Alfven 
speed, and consequently lowers the MRI threshold to below 
the Hall drift (e.g., for /s < 1). 

On the other hand, the circumplanetary disc may host 
a wide range of conditions. Here the Hall-MRI zone in the 
circumplanetary disc is resilient to changes in fdg, but re¬ 
cedes if Pz increases through enhanced density or reduce 
field strength, so that the field is predominantly toroidal. 
The field would be predominantly poloidal if Ohmic resis¬ 
tivity exceeds the Hall effect, either by enhanced resistivity 
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through higher density, or reduced Hall EMF from a lower 
held strength. However, further analysis into the potential 
for MRI turbulence in the circumplanetary disc will rely on 
models which are not affected by numerical viscosity. 


4 IMPLICATIONS FOR MHD GAP 
MODELLING 

In this section we consider the implications of including non¬ 
ideal effects in MHD modelling of gas how. We determine 
the extent to which diffusion limits small scale held gra¬ 
dients to calculate the minimum held bending radius. We 
discuss the implications of dihusion-limited held bending on 
the resolution needed for simulations. We also estimate the 
strength of large-scale magnetic forces, to compare with the 
hydrodynamical forces included in the TOM12 simulation. 
This allows us to determine where, if at all, magnetic forces 
have the most inhuence on gas how. 


Figure shows a contour plot of logjQ(Linin/au) at 
y = G, with logarithmically spaced contour level. The mini¬ 
mum held bending length-scale follows Lmin ex rjp/v every¬ 
where except the circumplanetary disc where the rjp is large 
and L„, L,, are small. It increases toward the planet, and 
is smallest at z ~ Hp in the protoplanetary disc where the 
dihusivity is lowest. 

This can be used to gauge the minimum resolution re¬ 
quired for MHD simulations. For example, as the minimum 
held bending length scale is large in the circumplanetary 
disc, the resolution is determined by velocity and pressure 
gradients length-scales instead. In the protoplanetary disc, 
signihcantly higher resolution is needed to resolve turbu¬ 
lence which can be tangled on a very small length scales. 
The Hall effect can pose an additional challenge to simula¬ 
tions as it is dissipationless and so can introduce small-scale 
held structure (e.g. Sano fc StonepOO^ Kunz Lesur|2013 


Lesur et al.||2014 O’Keeffe fc Downes||2014 Bai|2015 1. 


4.1 Minimum field-gradient length-scale 

Non-ideal ehects, particularly Ohmic resistivity and ambipo- 
lar dihusion, can wash out magnetic held gradients and limit 
the held bending length-scale, Lp- Here we calculate the 
smallest length-scale that the held can bend on, Lmin, given 
the level of magnetic resistivity and dihusion throughout the 
gap as calculated in Section [3.2| The role of the Hall ehect 
on magnetic held gradients is highly uncertain, being able 
to enhance or resist gradients. Therefore we neglect its ehect 
in this simple calculation, but note the potential key impact 
it can have on the large-scale structure we consider here. 


Just as equations (121, (141-( 151 specify the maximum 


dihusivity which permits the gradients needed for a tur¬ 
bulent or toroidal held, we can invert this relationship to 
specify the minimum bending length-scale given the mag¬ 
netic dihusivity. This is applicable to large-scale, quasi-static 
magnetic held structures (as opposed to rapidly huctuating 
turbulence), for which magnetic induction is balanced by 
magnetic dihusion. We derive the relation between magnetic 
dihusion and a general held geometry using the induction 


equation [equation (111]. While it is generally necessary to 
treat each component of the induction equation separately, 
we can estimate the dihusion limit for the minimum gradient 
length-scale. 


rjp < vLb 


L-^+Lfi 


(27) 


We calculate the velocity and dihusivity gradient length- 


i/ = //|V/|, 


(28) 


with / set to the relevant huid velocity or Pedersen dihu¬ 
sivity, respectively. Note that as we use the velocity in the 
frame orbiting with the planet to calculate L^, we remove 
the Keplerian component with respect to the star, which 


acts as a large constant ohset in equation (28l. 


We invert equation (271 to determine the minimum held 
gradient length-scale Lmin: 

-1 


Lmin — 


2rjp 

V 


1 - 


V 

vLp 


1 - 


JL 

vLr 


+ 


At) 

vLy 


(29) 


4.2 Relative strength of magnetic forces 

Large-scale magnetic helds have the potential to inhuence 
gas dynamics, rather than being merely passively drawn 
along with the gas. We touched on this in Section |2.4| 
when we compared the toroidal component of magnetic pres¬ 
sure with gas pressure. For example, circumplanetary jets 
have been launched in ideal and resistive MHD simulations 
( Machida et al.|2006 Gressel et al.|2013 ). 

The TOM12 simulation is not magnetised; nevertheless 
we can estimate the size of magnetic forces and compare 
them with hydrodynamic forces in the simulation. This will 
indicate where, if at all, they have the greatest potential to 
alter the balance of forces and so inhuence the gas how. It 
also acts as a consistency check, validating (or otherwise) our 
adoption of a hydrodynamic simulation as the underlying 
model. 

The full momentum equation, including both hydrody¬ 
namic and magnetic forces is 

1 dv o 

- — -I- (v • V) V = —2QpZ X V — QpZ X (z X r) — V'h 

--Vp -b ^ (V X B) X B, (30) 

p hvrp 

where r = {x, y, z) and is the combined gravitational po¬ 
tential from the star and planet: 


$ = - 


GM* 


GMp 


(31) 


The hydrodynamic terms in equation (301 are the time 
derivative and inertial force on the left hand side, and the 
Coriolis, centrifugal force, gravitational, and pressure gradi¬ 
ent forces on the right hand side. These are calculated in the 
frame co-orbiting with the protoplanet in the TOM12 sim¬ 
ulation. We neglect the time derivative because it is small 
with the simulation in almost steady state. 

The addition of the magnetic force [the Hnal term on 
the right hand side of equation @] will inevitably change 
the balance of forces, however the impact will be negligible 
unless the magnetic force is strong compared to the inertial 
force. We calculated the inertial force L/ = |(v • V) vj from 
the TOM12 simulation and its associated gradient length 
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X[au] X [au] x[au] 


Figure 4. Same as bottom-right panel in Fig.j^ except varying the dust-to-gas mass ratio f^g = 0,10“^ (top- and bottom-left panels), 
column density enhancement factor /s = 0.1,10 (top-and bottom-centre panels), and magnetic field enhancement factor fs =0.1,10 
(top- and bottom-right panels). 



-6 -4-2 0 2 


logioCLg /au) 
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Figure 5. Left panel shows contours of the minimum field bending-scale, Lmin? shown in the edge-on view at y = 0. Centre and right 
panels shows the strength of magnetic forces relative to hydro dynamical forces. The centre panel shows the logarithm of the ratio of the 
magnetic force to the inertial force, \ogiQ{FB/Fj) at y = 0. The dashed contour shows the location of Fq = Fj. The right panel shows 
vertical profiles of logiQ{FB/Fj) and \og-^Q{FB/Fp), the logarithm of the ratio of the magnetic force to the inertial and pressure gradient 
forces at a:: = 0.3 au, y = 0 as the solid and dashed curves, respectively. 
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scale 


L, = 


l(v-V)v|’ 


(32) 


using second-order finite differencing, with good convergence 
compared with first-order differencing. 

We estimate the magnetic force from large-scale non- 
turbulent fields using, 

F..^|(VxB)xB|^||. (33) 

We do this by nsing the Alfven speed calculated for the 
total, non-turbulent held strength. Our magnetic force esti¬ 
mate relies on a knowledge of the held gradient length-scale. 
Lb- We expect that, in general, Lb traces the huid gradient 
length-scales because either the held is dragged along with 
the huid or visa versa. If this is the case, we take Lb = Lx, 
and it cancels from the ratio of the magnetic and hydrody¬ 
namic forces Fb/Fi- If, instead, the gas varies over a shorter 
length-scale than dihusion permits the held to bend, we take 
Lb = Lmin as calculated in Section [4.1[ In this way, we com¬ 
pare the magnetic force over the same scale as the hydrody¬ 
namic forces, to the extent that diffusion permits. We have 
not included the the force arising from effective MRI turbu¬ 
lent viscosity, as it is already the subject of extensive study 
within the literature (see [Turner et aL]|2014[ ). 

We also consider the ehect magnetic forces have on the 
scale height by comparing them with the pressure gradient 
force, Fp = |(Vp)/p|, with length-scale Lp calculated using 


equation (28l. 


The centre-panel of Fig. shows a contour plot of the 
logarithm of the ratio of the magnetic force to the inertial 
force, logjQ(Fs/i^/) a.t y — 0. The dashed contour denotes 
Fi — Fb- The right panel shows the vertical proHle of the 
ratio of the magnetic force to the inertial, and pressure gra¬ 
dient forces. Magnetic forces are strongest in the disc atmo¬ 
sphere, exceeding the inertial and pressure gradient forces 
above 2 > 2.5-3.5Up, where they are able to influence the 
protoplanetary disc scale-height. 

Magnetic forces are weakest in the circumplanetary disc 
where strong magnetic diffusion limits the field bending 
length scale to the diffusion length scale. Although artificial 
diffusion lessens the reliability of the circnmplanetary disc 
structure in the TOM12 simulation, KW14 shows that mag¬ 
netic diffusion is strong across a range of circumplanetary 
disc models, owing to the column high density shielding the 
disc from external ionising radiation. The strong diffusion 
increases the minimum field bending length scale, thereby 
reducing the magnetic force. This indicates that large-scale 
magnetic forces cannot yield any significant accretion in the 
circumplanetary disc, which must come from the MRI or 


Gravitational Instability-MRI limit-cycles (e.g., Martin & 
Lubow|[2011 Lubow fc Martin|2012 1 . 


5 DISCUSSION 

In this paper we examined the importance of non-ideal ef¬ 
fects in determining the magnetic field structure in a gap 
snrrounding a giant protoplanet. We modelled the gap us¬ 
ing a snapshot from the pure-hydrodynamical gap-crossing 


use this snapshot as a basis to a posteriori estimate key MHD 
quantities semi-analytically, which would otherwise be very 
challenging to incorporate into simulations. We calculated 
the ionisation fraction produced by cosmic-rays, stellar X- 
rays and radioactive decay, including the effect of grains. We 
calculated Ohmic resistivity, ambipolar diffusion and Hall 
drift to determine whether an MRI field conld be generated 
and if a toroidal field could couple to the gas flow. We es¬ 
timated the magnetic field strength in the protoplanetary 
disc from inferred accretion rates, and determined the gap 
field strength from flux-freezing. 

We found that a magnetic field would be easily drawn 
from the protoplanetary disc into the gap and circumplane¬ 
tary disc. A toroidal field permeates the gap, but is weakened 
by expansion from magnetic pressure above a > 2-3 Hp. The 
gap is MRI unstable at 2 : > 0.5 F[p with tnrbulence extend¬ 
ing down to the midplane if the vertical component of the 
magnetic field is along the rotation axis (i.e., if s = 1). If 
however, the vertical field and rotation axis are anti-aligned, 
the conditionally unstable regions are non-turbnlent. As pro¬ 
toplanetary discs exhibit a range of field/rotation axis mis¬ 
alignment angles we expect the size of the MRI tnrbulent 
region to differ between systems, with significant implica¬ 
tions for the flow dynamics (Hull et al.|2013 Krumholz et al. 


20131. 


simulation by Tanigawa et al. (20121. Our approach was to 


This direction dependence of the MRI originates with 
the Hall effect. The Hall effect is strong below 2 < 2 Hp and 
plays an important role in generating MRI instability. We 
found that it can also facilitate conpling of a toroidal field, 
despite strong magnetic Pedersen diffusion, and influences 
the orientation of the toroidal component. This may lead to 
counter-wrapping of the toroidal field if s = 1, where the 
field drift opposes the Keplerian gas flow. Further simula¬ 
tions are needed to probe the role of the Hall effect in this 
system. 

By testing the sensitivity of the calculations to dust con¬ 
tent, column density, and magnetic field strength, we found 
that the gap was generally susceptible to the Hall- or ideal- 
MRI. On the other hand, the MRI in the circumplanetary 
disc may be quite sensitive to disc conditions and conld have 
a turbulent, toroidal or vertical field. As the circumplanetary 
disc evolves, the disc may experience a range of different field 
configurations in keeping with the varying column density 
and dnst content. Understanding the evolution of a circnm¬ 
planetary disc is key for satellite formation studies, which 
is believed to be the formation site of moons. Bimodality 
in circumplanetary disc dynamics caused by the Hall effect 
may transfer to the growth of moons within the system. For 
example, turbulent heating will effect the location of the 
circumplanetary ice-line. 

Finally, we have calcnlated the minimum magnetic field 
gradient length-scale limited by magnetic diffusion. Mag¬ 
netic diffusion resists field gradients, and so minimal field 
bending is permitted in the circnmplanetary disc where 
Ohmic resistivity is strong. We find that large-scale (non- 
turbulent) magnetic forces are unable to drive accretion in 
circumplanetary discs. Turbulence aside, we found the large 
scale-flow features are well modelled by a hydrodynamical 
fluid as large scale magnetic forces as small outside the 
protoplanetary disc atmosphere. Here, magnetic coupling is 
strong and so a strong field may be able to produce a jet, 
winds, or other variability. 
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In summary, we find that the magnetic field surround¬ 
ing a giant protoplanet is mostly toroidal, with large bands 
of ideal- and Hall-MRI turbulent zones. The magnetic field 
geometry is dependent on the orientation of the vertical field 
component, established during the collapse of the protostel- 
lar core. Non-ideal effects are important in the gap and need 
to be included in future MHD simulations of gap crossing. 
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APPENDIX A: WEAK-FIELD CONDITION 
FOR HALL-MRI 

Developing and sustaining MRI turbulence requires two con¬ 
ditions: (i) firstly, that field perturbations must grow fast 
enough, and (ii) secondly, the wavelength of the fastest- 
growing turbulent mode A, must be contained within a 
disc scale-height (i.e., A < H). Here we develop the weak- 
field condition for Hall-MRI. The wavelength of the fastest- 


growing Hall-MRI mode is [see equation (B14) of Wardle & 
Salmeron||2012| 


A= - 

V 


- 4,r]pv - + 




(Al) 


where v is the growth rate of the fastest-growing mode. In 
the Hall MHD limit, r]p =0, the maximum growth rat e 
attains the ideal rate v = (Wardle & Salmeron|2012 1. 

Using this result in strong-Hall limit, TyrrDuT,; 1, the 
wavelength of the fastest-growing mode is approximately 


\^2tt 


(A2) 


up to a constant factor of order unity. 

Ensuring that the two MRI conditions [(i), (ii) above] 
are met allows us to bracket the Hall-drift: 




^ Vh D 


(A3) 


which, with the thin-disc approximation H = Cs/tl, leads to 
the weak-field limit 


Pz = ^5— > Stt . 

Va,z 


(A4) 


This is identical to that in the ideal and resistive MRI 
regimes [equation (|13|); Okuzumi & Hiros^|201l|. 


APPENDIX B: NON-IDEAL EFFECTS ACTING 
ON A TOROIDAL FIELD 


Here we determine the coupling threshold for a toroidal 
field using the azimuthal component of the induction equa¬ 
tion. We consider a Keplerian disc with velocity, v = Qrcf), 
independent of height, and uniform transport coefficients 
'no,'nH,riA for simplicity. 

Shear generates a toroidal component from a poloidal 
field in the inductive term: [V x (v x B)]^ = — ^Q,Br. The 
resistive term, [—V x (rjoV x B)]^ = rio{V^B^ — B^/r'^) ~ 
rjoB^/, is dominated by vertical gradients (see Section 

, as are the Hall and amblpolar 
Comparing the two terms yields 

77 < (Bl) 


2.4| [Wardle & Konigl||l993 
terms, which scale similarly, 
the coupling criterior0 


where r] — rjp,rjH- 

How will the field evolve once it decouples from the gas 
motion? Recasting the induction equation to show the field 
line drift, Vb, makes this clear ( Wardle & Salmeron|20l'^ : 

aB 


dt 

where 


= V X |(v J- Vb) X B - 770 [(V X B) • b] b| , (B2) 


Vb = ^ (V X B) 


X B 


^(VxB)^ 


(B3) 


The azimuthal component of the drift velocity, Vb,^,, con¬ 
trols the wrapping of the field lines with the shear. Neglect¬ 


ing terms of OlHjr), the azimuthal drift velocity is (Braid- 
ing & Wardle|2012 1: 

77P BzBfj)^s Br^s 

ll~B~ 


= 


H R2 
= Vp + Vh, 


(B4) 

(B5) 


where we have separated the Pedersen and Hall drift com¬ 
ponents. Pedersen drift always resists winding, tending to 
reduce \B^^s\ and straighten field lines. On the other hand, 
the Hall component will enhance or resist winding according 
the sign of Br,s- 

In the natural configuration for s = 1 (i.e., Br,s > 0, 
and B^^s < 0 for z > 0), Vh < 0 and Vp < 0 so that the 
components cooperate in unwrapping the field. For domi¬ 
nant Hall, the unwrapping will overshoot, winding the field 
in the opposite direction so that B^^s > 0. For the anti¬ 
aligned field, s = — 1, the natural configuration is Br ,3 < 0, 
and > 0 above the midplane. Once again, Pedersen re¬ 
sists wrapping of the field lines and if it dominates the field 
will tend towards vertical. On the other hand, now Hall drift 
tends to enhance field motion along so that even if the 
field is somehow wound against the flow (i.e., B^^s > 0) both 
components cooperate to restore B,j ,^3 > 0 in equilibrium. 

This behaviour is summarised in Fig. |B1| which shows 


® This threshold is a factor of 10 (H/r) ^ ~ 5 X 10^ lower than 
[Turner & Sano] ( |2008| l, as we allow for vertical gradients in the 
non-ideal terms. 
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Figure Bl. Regions of r]jf-rjp parameter space, shaded accord¬ 
ing to the toroidal field drift. Ideal MHD (light green) and Hall 
dominated MHD (dark green) will sustain a toroidal field, al¬ 
though the orientation of the field varies with the vertical field 
orientation in the Hall-MHD region. Pedersen diffusion only per¬ 
mits a vertical field in the yellow region. 


how the resulting field geometry varies in the rjH-rjp pa¬ 
rameter space. Non-ideal effects are weak in the lower-left 
quadrant (light green), with the toroidal field wound up by 
the disc so that < 0. Pedersen diffusion dominates 

above the diagonal (yellow region), so that the field tends 
toward vertical. The Hall effect counteracts diffusion below 
the diagonal (dark green), so that the toroidal component is 
enhanced for s = — 1 or counter-wrapped when s = 1. 
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